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In this work, we systematically study the interaction of D* and nucleon, which is stimulated by the observation 
of A c (2940) + close to the threshold of D" p. Our numerical result obtained by the dynamical investigation 
indicates the existence of the D*N systems with J p = | , | , which not only provides valuable information to 
understand the underlying structure of A c (2940) + but also improves our knowledge of the interaction of D* and 
nucleon. Additionally, the bottom partners of the D'N systems are predicted, which might be as one of the tasks 
in LHCb experiment. 
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I. INTRODUCTION 

The BaBar Collaboration reported a new charmed hadron 
A c (2940) + with mass M = 2939.8 ± 1.3(stat) + l.O(syst) 
MeV/c 2 and width F = 17.5 + 5.2(stat) ± 5.9(syst) MeV 
by analyzing the D°p invariant mass spectrum, which is an 
isosinglet since there is no evidence of doubly charged part- 
ner in the D + p spectrum (H. Later, the Belle Collaboration 
confirmed A c (2940) + in £ c (2455) a++ 7r + '- channels Jl, which 
gave M = 2938.0 ± 1.3^.0 MeV/c 2 and T 
consistent with the BaBar's measurement U] 

The observation of A c (2940) + has stimulated extensive in- 
terest among different theoretical groups, which have pro- 
posed different explanations to the underlying structure of 
A c (2940) + . Since A f (2940) + is near the threshold of D*p, 
A c (2940) + is explained as an S-wave D*°p molecular state 
with spin parity J p = i , where the obtained decay behav- 
ior of A c (2940) + is not only consistent with the experimental 
measurement but also is applied to test the molecular struc- 
ture 0]. Later, the strong decay of A f (2940) + was studied 
under the D*N molecular state assignments with J p = 1 and 
J p = i + in Ref. H, which indicates that A c (2940) + should be 
assigned as the D*N molecular state with J p = \ + . Recently 
the radiative decay of A c (2940) + under the assignment of the 
D*N molecular state with J p — j + was performed in Ref. fH]. 



i + 



D*N molecular state, D* interacts 



If A c (2940) + is J' ■ 
with nucleon via P-wave. 

Besides giving these exotic explanations to A f (2940) + , the- 
orist has tried to find suitable assignment to A f (2940) + under 
the framework of the conventional charmed baryon. The po- 
tential model once predicts the masses of A c with J p = | , | 
are 2900 MeV and 2910 MeV H 01, respectively. Cheng 
and Chua calculated the ratio of H c nlY. c n if A c (2940) + is of 



J p = 2" 

J 2 



or | in heavy hadron chiral perturbation theory 



and indicated that such ratio is useful to distinguish the J p 
quantum number of A c (2940) + . In Ref. |9J], the strong decays 



of the newly observed charmed hadrons have been calculated 
by the 3 P model. The corresponding numerical result indi- 
cates that A c (2940) + could only be as D-wave charmed baryon 

ci (5 ) or ci ( 1-^ *- see t ' le notat i° ns f° r tne D-wave charmed 
baryons in Ref. J9Q) while A c (2940) + as the first radial excita- 



tion of A e (2286) + is fully excluded since A c (2940) 



D°p 

was observed by BaBar QlJ] . In the relativistic quark-diquark 
model, Ebert, Faustov and Galkin suggested A c (2940) + as 
the first radial excitation of Z c with J p — | + Holl . The re- 
sult obtained by chiral quark model indicates that A c (2940) + 
is D-wave charmed baryon A C 2 D^(| + ifllh . In Ref. Il2ll . 
A c (2940) + as the first radial excitation of the S c with J p = | + 
was proposed by solving the three-body problem by the Fad- 
deev method in momentum space. By the mass load flux tube 
model, the authors in Ref. 01] suggested that A c (2940) + 
could be as the orbitally excited A^ with J p — | . 
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FIG. 1: The observed decay modes of A c (2940) + and the comparison 
of the mass of A c (2940) + with the thresholds of D*p, Dp, 2^(2455)^. 
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Although different theoretical explanations to A c (2940) + 
were proposed, at present the properties of A f (2940) + are still 
unclear, which means that more theoretical efforts are needed 
to reveal its underlying structure of A e (2940) + . 

As shown in Fig. [1] A c (2940) + not only decays into Dp and 
£ c (2455)7r, but also is close to the threshold of D*p, i.e., about 
6 MeV mass different between A c (2940) + and the threshold 
of D*p. Thus, exotic D*N molecular state becomes one of 
the possible explanations to the structure of A f (2940) + . The 
dynamical study of the D*N system in one-boson exchange 
model is an interesting research topic at present, which can 
help us further clarify the D*N molecular state assignment to 
A c (2940) + and deeply understand the D*N interaction. In this 
work, we systematically carry out the dynamical investigation 
of the D*N system. 

This paper is organized as follows. After introduction, we 
present the detail of the dynamical study of the D*N system, 
which includes the relevant effective Lagrangian and coupling 
constants, the detailed derivation of the effective potential of 
D*N interaction, the corresponding numerical result, the study 
of B*N system. Finally, the paper ends with the discussion and 
conclusion. 



are two by two pseudoscalar and vector matrices 



V2 



V2 " V2 
P~ 



V3 V5 > 



By expanding Eqs. (Q]i-®> one further obtains the effec- 
tive Lagrangian of light pseudoscalar meson P coupling with 
heavy flavor mesons 



£** 



9": 



2g 

T 

2 8, 



ha- 



= -^(P b <p* a \ + <Pl A <pl)d A V ba . 

Jn 



(5) 
(6) 



The effective Lagrangian describing the coupling of light vec- 
tor meson V and heavy flavor mesons reads as 



£„y = -<2Bg v P b P\v' 



(7) 



-Cp"py 



)ba, 



Lp-rv = <2Bg v P\P* c J 



' ha 



+i2 ^AgvP^Vfid^v - dyVJba. 



(8) 



(9) 



II. THE DYNAMICAL STUDY OF D'N SYSTEM 



The effective Lagrangian of the scalar cr interacting with 
heavy flavor mesons can be further expressed as 



A. The effective Lagrangian and coupling constants 

In this work, we perform the dynamical study of D*N sys- 
tem. In order to deduce the effective potential of the D*N 
interaction resulted from the pseudoscalar, vector and scalar 
meson exchanges, we adopt effective Lagrangian approach. 
In this section, we collect the relevant effective Lagrangian. 

In terms of heavy quark limit and chiral symmetry, the La- 
grangians depicting the interactions of light pseudoscalar, vec- 
tor and scalar mesons with S-wave heavy flavor mesons were 
constructed in Refs. lfl4l-l20ll 



-Chhp 
-Chhv 



HHcr 



ig(H b7 ,Aly 5 H a ), 

+iA{H h cr llv F> lv (p)H a ), 
g s (H a crH a }, 



(1) 

(2) 
(3) 



where the multiplet field H is composed of pseudoscalar P 
and vector P* with P (t)T = (D w+ , D w °) or (Z? w °, And 
H is defined by 



[TV-Pars]. 



(4) 



Here, H = y #Vo and v = (1, 0). 

In the above expressions, the P and P* satisfy the normal- 
ization relations (0\P\Qq(0~)) = and {O^Qqil-)) = 
e^yfMp - . The axial current is A'' = \{g 'd t £-%d v g') = fd^+ 



■ ■ ■ with f = exp(iP//J and f„ = 132 MeV. (T ba = ig v Y M ba 

Fpy{p) = dtfy ~ dyPfi + fa, py], ATld g V = OTp/ f W - USK, I 



/V2, 
and 



Zrra = -2g s P b P\a, 



(10) 

(ii) 

As calculated in Ref. Ill, g = 0.59 in Eqs. Q-© is ob- 
tained by the full width of D* + determined by experiment. 
The parameter B appearing in the Lagrangians relevant to vec- 
tor meson can be fixed asB- 0.9 by vector meson dominance 
mechanism while A = 0.56 GeV 1 can be obtained by com- 
paring the form factor calculated by light cone sum rule with 
that obtained by lattice QCD. As the coupling constant related 
to scalar meson cr, g s = gj(2 V6) with g„ = 3.73 was given 
in Refs. HQ. 

The effective vertices depicting the interaction of nucleon 
with pseudoscalar meson P, vector meson V and scalar meson 
cr are respectively 



FNN 



gPNN 

^2m^ 



NbJsy^d^taNa, 



(12) 



■CvNN - 
-CirNN = 



- ^gvN N N b (y, + ^-cr^Y^Na, (13) 



gaNNNa(rN a , 



(14) 



where N T = (p, n) represents the nucleon field. The coupling 
constants g 2 „ NN l{An) = 13.6, g 2 pNN /(^) = 0.84, gl m l{An) = 

20 and gl m /(4n) = 5.69 with k = 6.1 (0) in Eq. fl3 for 
p (a>), which are used in the Bonn nucleon-nucleon potential 
rf23ll and meson productions in nucleon-nucelon collision 1I24I - 
12611 . We follow the convention for the signs of coupling con- 
stants as given in Refs. I24U26I1 . 

We need to emphasize that in this work we only consider n, 
p, ai and cr exchanges due to the weak coupli ng of tj or <p to 
nucleons as indicated in many previous works 1231 12411 . 
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B. Derivation of the effective potential of D'N interaction 

The scattering D'N — > D'N occurs via n, p/co and cr ex- 
changes. 

The scattering amplitude iM{ J, Jz), which is obtained by 
effective Lagrangian approach, is related to the interaction po- 
tential in the momentum space in terms of the Breit approxi- 
mation 



V(q) 



1 



vn,-2M,n/2M 



■M(JJz) 



(15) 



/ 



where M, and Mj denote the masses of the initial and final 
states, respectively. The potential in the coordinate space V(r) 
is obtained after Fourier transformation. For compensating 
the off-shell effect of exchanged particle and describing the 
inner structure of every interaction vertex, the form factor is 
introduced with monopole form F(q 2 ) = (A 2 - m?)/(A 2 - 
q 2 ) when writing out scattering amplitude, where the cutoff A 
should be around 1 GeV ll22tl . 

With the above preparation, the ji exchange potential be- 
tween heavy flavor meson D* and nucleon in the momentum 
space is obtained 



gnNN I 



-OF • q){<T ■ q)P(q 2 )F(q 2 ) t d ,-t n . 



m N 



(16) 



The p exchange potential can be written as 
1 



^>(q) = gpNNg 



e m ) 



1 + (1 + 2a-) 



icr ■ qxQ 



Am 2 



+A 



i—T ■ qxQ + — + —T x q ■ cr x q 



xP(q 2 )F{q-)T D , ■ tn. 



(17) 



The a> meson exchange potential can be easily obtained by 
replacing the relevant coupling constants and the mass of ex- 
changed light meson, and removing the isospin factor td- ■ 
and setting k — in Eq. dTTb . The cr exchange potential reads 
as 



VAq) 



gcrNNgsC 



1 



cr ■ qxQ 

Am 2 



Y(q 2 )F{q 2 ). 



(18) 



In the above expressions of the obtained potentials, T a = 
ie 0a P y ^p f e^' and P (q 2 ) = nl \ m2 ■ The polarization vectors are 



q-+m- 



defined as e* = + -±=(1, +;', 0), e° = (0,0, 1). Here, is the 
mass of the exchanged meson for D'N — » D'N transition. 



In this work, we focus on the D'N systems with the to- 
tal angular momentum / < |, which are of positive or neg- 
ative parity. Such D'N systems can be categorized as twelve 
groups according to the quantum number I(J P ) of system, i.e., 
the systems with O(^), 1(1*), 0(| ± ), Iff*), 0(f ± ) and 1(| ± ). 
Each of the D'N systems with J — \ is composed of two 
states 



I 1" \ I 2 \ I 4 \ 

1 i + \ 1 4 \ 1 

/(-)): P.), I 2 / 5 .). 



(19) 
(20) 



And each of the D'N systems with J — |, | is constructed by 
three states 



«|">): 








(21) 








fay- 


(22) 


«!*>>= 








(23) 


«!*>): 


H 


H 




(24) 



Here, we use notation 2s+l Lj to show the concrete informa- 
tion, which includes total spin S , angular momentum L, total 
angular momentum / of the D'N system. S , P, D, F and G in- 
dicate that the couplings between heavy flavor meson D* and 
nucleon occur via S -wave, P-wave, D-wave, P-wave and G- 
wave interactions respectively, which means that in this work 
we will include such /-wave contributions (i = S, P, D, F, G). 

The general expressions of these states in Eqs. (fT9b- (f24b can 
be explicitly written as 



|2S+1 



m,m' ,mi,nis 



(25) 



where C J , M , Cl M , and C Sms 



are Clebsch-Gordan co- 



m,Lmi Sm s ,Lm L i,„j m / 

efficients. Yi„, L is spherical harmonics function. xi m denotes 
spin wave function of the corresponding state. The polar- 
ization vector for D* is defined as e™ = +-^(f'" ± ity) and 
e ( "' = el"- Here, the polarization vector in Eq. dZ5b is just the 
one appearing in the potentials listed in Eqs. ([T6l)-(IT8"1). 

According to the sub-potentials in Eqs. (IToTl-lfTSl) and wave 
functions in Eqs. (fT9b-<f25b. one obtains the total potential of 
the D'N system with J p = \ ± , 
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Vi- 



-D-2C 
y/2T 



-D + C 



yf2T 
180' 



0-2T 



-D-2C + 80' + \0 y/2T + 2 V20' + ^0 
yf2T + 2 V2C + ^fo -D + C+ 100' - f O - 2T 



(26) 
(27) 



which are two by two matrixes since the D*N system with 



or F 



1 + 



is constructed by two states just listed in 



Eq. ( fT9b or (f20l >. Analogously, the total potential for the D*N 



system with J p 
matrix as 



4 or ^ 



can be expressed by three by three 



V 3 - = 



c 



2T 



-T 



2T C-§£>- fO + 120' §D + jfO + 60' + r 



-T \D 



\0 + W + T -D-2C+ 120' 



15°) 

(-6D + 420' +40- 3T) 
^2 (-6D + 420' +40-37) C - ^(-12D + 80 + 210' + 15T) 



-2C-D-W + |0 



-2C -D-40' + \0 -^(-60 + 30' + 20 - 3r) 

3 3 V5 v ' 

^(-6D + 30' + 20 - 3F) C - ^(-3£> + + 150' + 6T) 



2^T 

4 V6 t^ 
7 J 

C + 4 [-9D + 5(840' - 30 - 4r)] 



C - 6£> - 40' + |r 



-Vf 3 

4 V6t 



C - ±(3D - 900' + 40 + 6T) ) 

4%/6 x 



-2C-D + 160'- ±0 



^=[15D + 7(300' + 50 + 3T)] C - j^D + 140' - §0 



^[150 + 7(300' + 50 + 3D] 

25 n , 85 /O _i_ 1()T 



(28) 



(29) 
(30) 

(31) 



with 



D = Dfrjy ■ t n + Dt - D' a , 

C = -C„- 2[(1 + k)C a t d , ■ t n + Cil 

T = -T„ + 2[(l + k)T a t d - ■ t n + T A ], 

O = a -[{\+2K)O S p T D >-T N + O 3 J, 



O' 



A p T D . -T N + A 



(32) 
(33) 
(34) 
(35) 
(36) 



where the expressions of C,, Z)', T t and O, (; = n, cr) are 
defined as 



Ti = 



Z (mi, r)Y Q (mi, r) - Z (A, r)Y (A, r) 



1 

— 

4n 

-|-(Ar-2)F (A,r) 



(38) 



r (m,-, r)Fo(mi, r) - 7q(A, r)F (A, r) 



1 

4^ 
-|-(Ar+l)F (A,r) 



(39) 



O, = 



47T 



O (m ; , r)Y (m h r) - O (A, r)F (A, r) 



-^-rF (A,r) 



(40) 



O' = 4mlA 



4m' 



F (mi,r)-Fo(A,r) 



(37) 



The expressions of C 1 ., D. ] , Tj and 0\ {i = p, co and j = ft, A) 
are similar to those of C,-, D'., Tj and O, respectively, which 

can be easily obtained by replacing factor Fj with Fj in Eqs. 
03)-(|40l) (F„ = gg„ m /(6 j2m N f„l Fy = gvNNgvPl^m 2 N ), 
Fy = gvAwg^/(3mjv), F a = go-NNgo-/(4m 2 N )). Here, 
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= ^A 2 - m 2 and functions Yq(x, r), Zq(x, r), Tq(x, r) and 
Oq(x, r) are 



solving the coupled-channel Schrodinger equation. One uses 
the FESSDE program fl27ll to produce the numerical results. 



Y Q (x, r) = , Z (x, r) - x 

r 

T (x,r) = j(2[l + l + _L_l 
L xr (xr) L 1 

O (x,r) = x 2 [l + -i-l. 

L xr (xr) 1 J 



C. Numerical result 

As shown in Se c. JHI the positive values of parameters g, 
PIAmdg s EE 

|21l|22|], which are relevant to pseudoscalar, 
vector and scalar meson exchanges respectively, are widely 
adopted in previous theoretical work. Thus, we first present 
the numerical result under taking the positive values of pa- 
rameters g, pi A and g s . 

With the isovector D*N system as an example, one shows 
the line shapes of C„ D„ T h <9„ C\, D j ., Tj, 0\ in Fig.|2]under 
taking cutoff A = 1 GeV, where the signs of g, pi A, g s are 
taken as positive. The pion exchange really provides impor- 
tant contribution to the exchange potential as shown in the left 
figure while the vector meson exchange also give consider- 
able contributions to the effective potential. The contribution 
from <x> meson exchange is comparable with that from pion ex- 
change due to the large coupling constant relevant to oj. Com- 
pared with 7i, p and oj exchanges, scalar meson exchange only 
gives small contribution to the effective potential. The spin- 
orbit terms O, from the vector and scalar meson exchanges, 
which correspond to the relativistic correction, are neglectable 
compared with the other terms. 
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FIG. 2: (Color online). The dependence of C%, A, Tu Oi, C\, D^, 
Tf, 0\ on r for the D*N system. Here, we use the D*N system with 
isospin / = 1 as an example. The cutoff A = 1 GeV and the signs of 
g, PI A, g s are taken as positive. 



Using the potential obtained above, the binding energy for 
"!>§>§ can be obtained by 
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FIG 3: (Color online.) The A dependence of the binding energy, 
the root-mean-square radius r and P\ which is the possibility of the 
\ 2 S i >, | 4 />i >, | 4 5 3 > and | 2 P3> states in the D*N systems with ^, | + , 

rand. 34 " 



2 , respectively. The thick dashed, the thick dash-dotted, the 
solid and the dashed lines are the bound state solutions of the D'N 
systems with I(J P ) = 0(f ~),0(| + ), 1(|~), 1(| + ) respectively. 



The obtained binding energy and the relevant root-mean- 
square radius r (in the unit of fm) of the D*N systems are 
presented in Fig. [3] with the variation of the cutoff A in the 
region of 0.8 < A < 1.2 GeV. Here, we only show the 
bound state solution with binding energy less than 10 MeV 
since the OBE model is valid to deal with the loosely bound 
hadronic molecular system. As shown in Fig. [3] we can 
find bound state solutions only for four D*N systems with 
I(J P ) = 0(|"), 0(^ + ), K5 ), 1(| + ) among twelve systems 
shown in Eqs. dT9b-324t. In Fig.[3j one also presents the possi- 
bilities of \ 2 S 1 ), \ 4 P 1 >, \ 4 S 3 ) and \ 2 Pi ) components appearing 

in the corresponding D*N systems with i , | + , | and | + , 
which indicates that \ 2 S 1), | 4 Pi), | 4 5i) and \ 2 Pi) states are 

2 2 2 2 

3- 
2 



and I respec- 



the D*N systems with J p = 



dominant in the D*N systems with j , j 
tively. 

As an isoscalar state, A ( l :(2940) can directly correspond to 
one of the D*N systems with I(J P ) = Q(^), 0(| ± ), 0(f *). 
If A+(2940) is D*N molecular state, such D*N molecular 
state should be of the binding energy -6 MeV. Among the 
above six possible D*N systems, only two D*N systems with 
I(J P ) = 0(i + ), 0(| ) are of -6 MeV binding energy under 
taking cutoff A as 1.17 GeV and 1.09 GeV, respectively. Thus, 
our dynamics study presented in this work support to explain 
A+(2940) as D*N system with I(J P ) = 0(| + ) or 0(|"). 

Although positive values for parameters g, pi A, g s are 
adopted in former work Il8l Il9l l2lL 12211 corresponding to the 
case with " + + + " in Fig. |H in fact the signs of these pa- 
rameters can not be well constrained by the experiment data 
or theoretical calculation, which could results in changing the 
signs of corresponding pion, vector and sigma exchange po- 
tentials of D*N systems. As shown in Fig. [4] one presents the 
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FIG. 4: (Color online.) The binding energy for the D'N systems. 
The solid line, the dashed line, the dash-dotted line, and the dotted 
lines are for the bound state solutions with J p = ~ , | , ~ , | re- 
spectively. Here, +/- in " + 1 ± 1 ± 1" denotes that we need to 
multiply corresponding pion, vector and sigma exchange potentials 
of the D'N systems listed in Eqs. dl6|>-l|18l> by an extra factor + 1/ - 1, 
which come from the changes of the signs of coupling constants. 
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FIG. 5: (Color online.) The A dependence of the binding energy, 
the root-mean-square radius r and P\ which is the possibility of the 
\ 2 S i >, | 4 Pi >, | 4 S 3 > and \ 2 Pi) states in the B*N systems with ±~, | + , 
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solid and the dashed lines are the bound state solutions of the B*N 
systems with I(J P ) = 0(f ~),0(| + ), 1(±~), 1(| + ), l(f + ) respectively. 



binding energy dependent on A under eight combinations of 
the signs of g, pi A, g s , where +/- denotes that we need to 
multiply corresponding pion, vector and sigma exchange po- 
tentials of D*N systems listed in Eqs. (fT6b-(n~8T> by an extra 
factor +1/ - 1 which results from the changes of the signs of 
8,0/ 1 8m- 

We find that the sigma exchange contribution can be negli- 
gible since the line shapes of binding energy dependent on A 
shown in the second column almost keep the same as those in 
the first column as describing in Fig. [4] n and p/ai meson ex- 
changes play very important role to form the D* N bound state. 
These observations are consistent with the expected behavior 
of the potential of D*N. 

Our numerical results shown in Fig.|4]indicate that there do 
not exist D*N molecular states with J F = | . 

D. Bottom partner of the D'N system 

In the previous section, the D*N system are investigated and 
the bound states solutions are found. The observed A c (2940) + 
can be assigned as the D*N molecular state with 0(4 ) or 

0(| ) supported by the dynamics study of the D*N system. 
Due to the heavy quark symmetry, we also extend the same 
formulism to the B*N system, which is the bottom partner of 
the the D'N system. The numerical results for the B*N system 
are listed in Figs.[5]and[6] which are similar to the discussion 
for the D*N system. 

If A c (2940) + as D*N molecular state with I(J P ) = 0(^ + ) 
or 0(| ), one expects that there should exist the correspond- 
ing bottom partner of A c (2940) + . As indicated in Fig. |5j 
one indeed finds the bound state solutions of B*N systems 
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FIG. 6: The binding energy for the B* N systems. The solid line, 
the dashed line, the dash-dotted line, and the dotted lines are for the 
bound state solutions with J p = | , | ,| , | respectively. Here, 
+ /- in (±1 ± 1 ± 1) denotes that we need to multiply corresponding 
pion, vector and sigma exchange potentials of B*N systems listed 
in Eqs. dl6ll-lll8ll by an extra factor +1/ - 1, which come from the 
changes of the signs of coupling constants. 

i + ^ — — 

with 0(j ) and 0(4 )■ Thus, the experimental search for B*N 

molecular state will be helpful to deep our understanding of 

the underlying structure of A f (2940) + , which might be as the 

task in LHCb. 



III. CONCLUSION AND DISCUSSION 

Stimulated by the observation of A f (2940) + US], which 
is close to the threshold of D*p, we study the interaction of 
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D* meson with nucleon N, where the OBE model is applied 
to obtain the effective potential of D*N system. By solving 
Schrodinger equation, one can find the bound state solution 
of the D*N system, which will be helpful to answer whether 
there exists the D*N bound state corresponding to A c (2940) + . 

As indicated by the obtained numerical result, there exist 
bound state solutions (-6 MeV binding energy) for the D*N 
system with I(J P ) = 0(| ), 0(| ) taking reasonable cutoff 
A, which indicate that it is possible to explain A c (2940) + 
as isoscalar S-wave or isoscalar P-wave D*N molecular by 
performing the dynamical study. Additionally, we find the 
bound state solutions for the isovector D*N systems. Search- 
ing isovector D*N states might be as the task in future experi- 
ment. In this work, we also predicted the existence of the bot- 
tom partners of the D*N systems as the extension of the study 
of the D*N system. Carrying out the experimental search for 
the B*N molecular states will be an interesting topic, espe- 
cially for LHCb experiment. 

Searching for exotic nuclei is a very important research 
topic in hadron physics and nuclear physics, which not only 
helps us to understand the interaction of meson or hyperon 
with nucleon but also provides the important information to 
reveal some underlying problems in astrophysics. There are 
extensive studies of hypernucleus |28 31|. 77-mesic nucleus 
i32ll33tl . It is natural to expect the exotic nuclei composed of 
vector heavy flavor meson and nucleon. Experimental search 
for exotic nucleus composed of a vector heavy flavor meson 
(Qq or Qq meson with Q = c or b) and the nucleon might be 



as the main task at J-PARC, RHIC and FAIR |0 HI since 
the heavy flavor meson {Qq or Qq meson with Q = c or b) 
can be produced in a nucleon-rich environment, which also 
provides another approach different from the production pro- 
cess of A f (2940) + to test the our prediction of the D*N/B*N 
molecular states to some extent. 
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